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The suppression pattern of quarkonia in AA(nucleus-nucleus) collisions (where the 
plasma of deconfined quarks and gluons, the Quark-Gluon Plasma (QGP) formation is 
expected), along with the comparative quarkonia production in pp collisions provides im- 
portant understanding into the properties of the produced medium. Experimental results 
of Relativistic Heavy Ion Collider(RHIC) indicate a suppression on charmonium produc- 
tion in AA collisions. Muons from the decay of charmonium resonances are detected in 
the ALICE Experiment at the Large Hadron Collider(LHC) for pp and Pb-Pb collisions 
with a muon spectrometer, covering the forward rapidity region 2.5< y <4.0. Analysis 
of the nuclear modification factor (Raa) at forward rapidity are presented and compared 
with mid-rapidity results from electrons in the central barrel covering \y\ <0.9. The roles 
of suppression and regeneration mechanisms are discussed, as well as the importance of 
the results of the forthcoming p-Pb data taking for the estimate of cold nuclear mat- 
ter effects on quarkonia. Perspectives for the bottomonia measurements are also given. 
Quarkonia results via muon channel from CMS experiment at LHC are compared with 
ALICE quarkonia measurements. 



I. INTRODUCTION 

Quantum Chromo-Dynamics(QCD) is the 
well-established theory of strong interactions. 
Quarks are the basic constituents of QCD, 
which interact through the exchange of glu- 
ons. A thermalized system where the proper- 
ties of the system are governed by the quarks 
and gluons degrees of freedom is called the 
Quark — Gluon Plasma(QGP) Lattice 
QCD predicts this new phase, QGP, at high 
temperatures(~155-160 MeV) The mo- 

tivation of the relativistic heavy ion physics is 
the experimental study of the hadronic matter 
under extreme conditions of temperature pj}. 
This is the main reason for the present 
heavy ion program ongoing at the Relativis- 
tic Heavy- Ion Collider (RHIC, at Brookhaven 
National Laboratory in New York) and at 
the Large Hadron Collider (LHC, at the Eu- 
ropean Organization for Nuclear Research in 
Geneva). 

A wealth of ideas have been proposed in the 
past few decades on the experimental and the- 



oretical understanding of QGP properties 0- 
Hj]. The plasma is short lived and all sig- 
nals from QGP have background from the hot 
hadronic phase that follows the QCD phase 
transition. Therefore, it is important to un- 
derstand the physics in hot hadronic phase 
with different experimental probes Q. 



II. QUARKONIA IN HEAVY-ION 
COLLISIONS 



Among the possible probes of the QGP, 
the interest in heavy quarks is motivated by 
their unique role in the diagnostics of the 
highly excited medium created in relativis- 
tic heavy-ion collisions. Both experimentally 
and theoretically for over two decades, the 
properties of heavy quarkonium states (which 
are bound states of heavy quark-antiquark 
pairs, charmonium and bottomonium) in a hot 
and dense QCD medium have been intensely 
studied [H, |l(J. Heavy-quark mass is large 
and so heavy-quark production is believed to 



be occur largely within the earliest phase of 
the collision. Therefore, the measurement of 
quarkonia is expected to provide essential in- 
formation on the properties of the strongly- 
interacting system formed in the early stages 
of heavy- ion collisions . 

If a J /if) particle(a cc bound state) is placed 
in QGP, the colour charge of the charm quark 
c will get screened by the quarks, anti-quarks 
and the gluons of the plasma. The char- 
monium production was discussed as an im- 
portant probe of plasma [H, EH , along with 
the suggestion by Matsui and Satz 14| which 
affirmed that with increasing centrality in 
heavy-ion collisions a suppression of the J/ ip is 
expected. The dissociation of the charmonium 
bound-state in a dense medium is connected 
with the Debye screening of the binding po- 
tential. This is influenced by a deconfinement 
of color charges and thus intimately connected 
to the formation of a QGP [15l42l| . 

Comprehensive experimental results at 
SPS [jj [22] (including feeddown from other 
less bound resonances like ip(2S) and Xc) 
and RHIC of J/ip production in AA colli- 
sions clearly indicate that even the stro ngly 
bound J ftp ground state is suppressed (23l - [26| . 
However at LHC energies the J / ip production 
could be even enhanced due to the coalescence 
of uncorrelated cc pairs in the medium [§3|. 
Initial state effects like modifications of the 
parton distribution functions in the nucleus 
relative to the nucleon shadowing need to be 
taken into account. The final state effects 
like the nuclear absorption are expected to be 
practically insignificant at the LHC ener gies . 
Various statistical and transport models [271 — 
HO] are proposed for LHC energies. Studying 
the pA collisions at the LHC energies is cru- 
cial to quantify the role of initial shadowing 
effects. 



A. J/tj) R AA 

A Large Ion Collider Experiment (AL- 
ICE) [3l| is a general-purpose heavy-ion ex- 
periment to study the nature of the quark 
matter under extreme temperature (> 0.3 
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FIG. 1: (Color Online) Opposite sign di-muon in- 
variant mass distribution for the 0%-90% most 
central collisions for 2.5< y <4, integrated over 

PT- 



GeV) [32, |33| created in nucleus-nucleus col- 
lisions at the LHC. In the framework of the 
ALICE physics program, the goal of the Muon 
spectrometer is the study of quarkonia pro- 
duction [34-36], along with open heavy fla- 
vor production and low mass vector meson 
properties [37j via the di-muon decay chan- 
nel in pp, pA and A A collisions. Using the 
particle identification potential of the cen- 
tral barrel detectors (\rj\ <0.9) J/tp have been 
also measured in di-electron channel at mid- 
rapidity (Mil. 

The forward rapidity(2.5< y <4) quarko- 
nia analysis with dimuons in Pb-Pb collisions 
at -y/sNN = 2.76 TeV is based on a sam- 
ple corresponding to an integrated luminos- 
ity of 70 /ib -1 which corresponds to 17.7M 
triggered events. The invariant mass spec- 
trum of /J, + fi~ candidates for 0%-90% most 
central events is shown in Fig. [TJ The signal 
is extracted by performing a combined fit for 
the background and signal contribution. The 
background is described by a Gaussian with 
a width(cr) which varies as a function of the 
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FIG. 2: (Color Online) Inclusive J/ip Raa for 
mid-(blue) and forward(red) rapidity as a function 
of the number of participating nucleons measured 
in Pb-Pb collisions at = 2.76 TeV. 



FIG. 4: (Color online) J/ip R AA for high- 
Pt(5<Pt<8 GeV/c) in Pb-Pb collisions at ^/snn 
= 2.76 TeV as a function of the number of par- 
ticipating nucleons and compared with CMS re- 
sults [H. 
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FIG. 3: (Color Online) J/V Raa in Pb-Pb col- 
lisions at ^/snn = 2.76 TeV recorded in 2011 as 
a function of the number of participating nucle- 
ons and compared with PHENIX results [3S| in 
Au+Au collisions at Vsnn = 200 GeV. 
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FIG. 5: (Color online) Inclusive J /if) Raa as a 
function of the J/ip pr for 2.5<j/<4.0 and cen- 
trality 0%-90% compared with CMS !39] and 
PHENIX El data. 



mass value. For the signal description a mod- 
ified Crystal Ball function was used which is 
a convolution of a Gaussian and power law 
functions that can fit the tails of the measured 
signal. 

The nuclear modification factor(i?AA) al- 
low us to quantify the medium effects on the 
production of J/4 1 - Raa gives the deviation 
in J/ip yield in AA collisions relative to the 



scaled (according to the number of binary 
nucleon-nucleon collisions) yield of J / ip in pp 
collisions. The raw yield has been corrected 
for acceptance and efficiency. 

Figure [5] shows the inclusive J / ip Raa for 
mid- and forward rapidity for Pb-Pb colli- 
sions at ^/snn = 2.76 TeV, as a function of 
the average number of nucleons participating 
to the collision (iVp art ) which has been calcu- 
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FIG. 6: (Color online) Comparison of Raa as a 
function of the J/ip pr for 2.5<y<4.0 in two cen- 
trality classes of 0%-20% and 40%-90%. 
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FIG. 7: (Color Online) Inclusive J/ip Raa as 
a function of the number of participating nucle- 
ons measured in Pb-Pb collisions at ^/snn = 2.76 
TeV. These results are compared with statisti- 
cal hadronization model 28] along with transport 
models by Zhao et al. [29| and Liu el al. [30| re- 
spectively. 



lated using the Glauber model [40j . Its impor- 
tant to note that in ALICE for both the di- 
muon and di-electron channel the J/tp produc- 
tion can be measured down to px — GeV/c. 
The di-muon data sample analysed and shown 
here are for nine centrality classes from 0%- 
10%(central collisions) to 80%-90% (peripheral 
collisions). For the dielectron analysis the 



0%-90% 
AA 



0.497 ± 0.006(stat.) ± 0.078(syst.) 



indicate J/ip suppression, improving the sta- 
tistical significance of the Raa measurements 
obtained with 2010 data |36|. The centrality 
integrated mid-rapidity value is i?^ _80% = 
0.66±0.10(stat.)±0.24(syst.). The systematic 
error is dominated by the pp reference both for 
mid- and forward rapidity measurements. 

Comparison with the RHIC measurements 
at = 200 GeV [H from the PHENIX 

experiment shown in Fig. [3] exhibit that the 
inclusive J/ip R A a at 2.76 TeV in the AL- 
ICE forward rapidity region are higher than 
that measured at 200 GeV in the rapidity do- 
main of 1.2 < \y\ < 2.2. High-pT J/ip Raa 
measured in ALICE for 5<pt<8 GeV/c range 
is compared with CMS measurements at cen- 
tral rapidity(|y| <0.9) for 6.5<p T <30 GeV/c 
in Fig. |U For ALICE the J/ip Raa at most 
central collisions is ~ 0.35 and the integrated 



value of R 



0%-90% 
AA 



0.384 ± 0.014(stat.) ± 



0.074(syst.). We observe from Fig. [4] that se- 
lecting high-pT drives down the Raa ■ 

Figure [5] shows Raa as a function of pt 
as measured by ALICE(red) results are com- 
pared to the CMS data [39l| (bl ue) as well as 
results from the PHENIX [38| (black) at a 
lower collision energy. The Raa is decreas- 
ing from 0.6 for low pr to about 0.4 at higher 
p T . The CMS results (0%-100% centrality, 
1.6 < \y\ < 2.4, pt > 6.5 GeV/c) are in agree- 
ment with the ALICE measurements (0%-90% 
centrality, 2.5<y<4, p? > 0) in the overlap- 
ping transverse momentum range whereas the 
lower energy results from PHENIX (0-20% 
centrality, 1.2 < \y\ < 2.2) show a signifi- 
cantly smaller Raa- Figure [5] shows the pr 
dependence of Raa obtained in the most cen- 
tral (0%-20%) and most peripheral (40%-90%) 
centrality classes. The J/ip suppression pat- 
tern depends on the centrality of the collisions. 
In the most central collisions (0%-20%), the 
suppression increases with the transverse mo- 
mentum of the J/ip. For the peripheral col- 
lisions (40%-90%), the pt dependence of the 
Raa is weaker and compatible with a flat be- 



haviour. Thus, for pt > 3 GeV/c the J/tp 
suppression is larger in central collisions. 

At forward rapidity, the non-prompt J/tp 
was measured by the LHCb collaboration to 
be about 10% in pp collisions at \/s= 7 
TeV [41] in the px range of the present analy- 
sis. Neglecting the shadowing effects and con- 
sidering the scaling of beauty production with 
the number of binary nucleon-nucleon colli- 
sions, the prompt J/ip Raa is estimated to 
be(upper limit), 11% smaller than the inclu- 
sive measurement. While estimating the in- 
fluence of non-prompt J/tp as a function of 
Pt and y on the inclusive Raa results, the 
LHCb measurement at y/s — 2.76 TeV is in- 
terpolated using CDF and CMS data. J/tp 
from beauty hadrons have a negligible influ- 
ence on the present measurements while as- 
suming a range of energy loss for the b-quarks 
from .RaaO) = 0.2 to #aaO) = 1. 

Figure [7] compares the ALICE results 
with results from a statistical hadronization 
model [28l . as well as two different transport 
models |29|, H3] • The models describe the data 
within uncertainties for -/V part larger than 70. 
The Statistical Hadronization Model [28| as- 
sumes deconfinement and a thermal equili- 
bration of the bulk of the cc pairs. Then 
the charmonium production occurs at the 
phase boundary by statistical hadronization 
of charm quarks. The two transport model 
from Zhao et al. [29| and Liu et al. [3(| dif- 
fer mostly in the rate equation controlling the 
J/ip dissociation and regeneration. Both the 
transport model calculations are shown as a 
band which connects the results obtained with 
(lower limit) and without (higher limit) shad- 
owing, which can be interpreted as the un- 
certainty of the prediction. The model from 
Zhao et al. incorporates a simple shadowing 
estimate leading to a 30% suppression for the 
most central Pb-Pb collisions assuming the 
charm cross-section da cS /dy « 0.5 mb at y 
— 3.25. The J/tp from beauty hadrons is es- 
timated at 10% and no quenching is assumed. 
The model from Liu et al. has shadowing 
from EKS98 and a smaller charm cross-section 
da c5 /dy ~ 0.38 mb is used. 
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FIG. 8: (Color Online) Opposite sign dimuon in- 
variant mass distribution for the 40%-60% most 
central collisions and 0<pt<3 GeV/c range in 
2.5< y <4. 
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FIG. 9: (Color Online) ALICE results for 
(N^CiSj/Nj/^PbPb/iN^^/Nj/^pp obtained in 
0<pt<3 and 3<p T <8 GeV/c for Pb-Pb collisions 
at =2.76 TeV recorded in 2011. ALICE 

results are compared with CMS [13] values for 
3<Cp T <30 GeV/c and 1.6< y <2.4. 



B. tP(2S) 

The study of tp(2S) provides an interesting 
comparison with J/tp production [42I l43j . In 
particular, a substantial fraction of the J ftp's 
is known to originate from tp(2S) and Xc de- 
cays 44]. Examining the ratio of ip{2S) to 
J/tp production can render a further insight 
on the mechanism affecting quarkonium in a 
hot and dense medium which was studied in 
detail at SPS energies [45|,|46j]. In a scenario, 
in which quarkonium states are suppressed 
by a Debye screening mechanism, the tp(2S) 
meson melts at lower temperatures, being a 
less bound state with respect to the J/tp }42| . 
A reduction of the tp(2S) nuclear modifica- 
tion factor, with respect to the J/tp can infer 
on this sequential melting (42J. However this 
picture might be complicated by charmonium 
production via recombination mechanisms at 
LHC, which can affect both the J/tp and the 
tp(2S). Production via recombination mecha- 
nisms should contribute mostly in the low-px 
region which can be studied by the ALICE 
experiment. 

The invariant mass spectrum of (J. + fJ,~ 
candidates for 40%-60% most central 
events and 0<pt<3 GeV/c is shown in 
Fig. [5] Figure M shows the double ratio 

( N i J (2S)/ N J/^)pbPb/(N 1 p { 2S)/ N J/i J )pp, which 
compares the ratio of tp{2S) over J/tp yields 
in Pb-Pb and pp. The pp results used here 
are from y/s = 7 TeV data-sample due to 
the low integrated luminosity at y/s = 2.76 
TeV. ALICE results which do not show 
large enhancement in central collisions are 
obtained in 0<p*r<3 GeV/c and 3<pt<& 
GeV/c for Pb-Pb collisions at y/sj^ = 2.76 
TeV and compared with CMS [47| values 
for 3<p T <30 GeV/c and 1.6<y<2.4 as a 
function of centrality. SPS measurements 
have indicated that tp(2S) is more suppressed 
than J/tp [Hj]. 



III. SUMMARY AND OUTLOOK 

Nuclear modification factor Raa of J/tp 
have been measured in Pb-Pb collisions at 



y/s^m — 2.76 TeV at forward rapidity and 
compared with mid-rapidity. It was shown 
that Raa decreases with increasing px and 
towards larger rapidity. Several models that 
take into account a recombination of cc-pairs 
are compared with the ALICE results. 

The inclusive J/tp production in Pb-Pb col- 
lisions at y/ snn = 2.76 TeV show less suppres- 
sion compared to PHENIX results at RHIC 
energies. The regeneration mechanism explain 
these experimental results. But a more quan- 
titative calculation requires the study of cold 
nuclear matter. There the roles of the sup- 
pression and regeneration mechanisms will be 
disentangled by quantifying the initial shad- 
owing effect with data from the forthcoming 
p-Pb run at LHC, planned in early 2013. 

Studying the bottomonia (T(15), T(25), 
T(35)) states and their suppression are impor- 
tant at LHC energies. Since the T is smaller 
and has larger binding energy than the J/tp, 
the states will provide valuable information 
complementary to the study of charmonium. 
The smaller bb rate results in a lower proba- 
bility for production by coalescence. In terms 
of A A collisions, the T is expected to disso- 
ciate at a higher temperature [48| - |50| than all 
the other quarkonium states, thus proving to 
be a more effective thermometer of the sys- 
tem mm- 
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